A quantum gate array can be programmed to evaluate the expectation value of any 

operator 
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A programmable gate array is a circuit whose action is controlled by input data. In this letter 
we describe a special-purpose quantum circuit that can be programmed to evaluate the expectation 
value of any operator O acting on a space of states of N dimensions. The circuit has a program 
register whose state \9(0))p encodes the operator O whose expectation value is to be evaluated. 
The method requires knowledge of the expansion of O in a basis of the space of operators. We discuss 
some applications of this circuit and its relation to known instances of quantum state tomography. 



An important feature of classical computers is that 
they can be programmed. That is to say, a fixed uni- 
versal device can perform different tasks depending on 
the state of some input registers. These registers define 
the program the device is executing. Quantum comput- 
ers II have a rather different property: Thus, Nielsen 
and Chuang established in [2( that a general purpose pro- 
grammable quantum computer does not exist. Such a de- 
vice would have to have the following features: It should 
consist of a fixed gate array with a data register and a 
program register. The array should work in such a way 
that the state of the program register encodes the unitary 
operator U that is applied to the state of the data register. 
As shown in 0], such devices cannot be universal since 
different unitary operators require orthogonal states of 
the program register. However, some interesting exam- 
ples of programmable devices could still be constructed. 
For example, non-deterministic programmable gate ar- 
rays were first considered in 2\ and analyzed later in a 
variety of examples [j| . More recently, quantum "multi- 
meters" were introduced and discussed in Jj]. Such de- 
vices are fixed gate arrays acting on a data register and 
a program register, together with a final fixed projective 
measurement on the composite system. They are pro- 
grammable quantum measurement devices |4| that act 
either non deterministically or in an approximate way 
(see 0). 

In this letter we will describe a different kind of pro- 
grammable quantum gate array that is useful to solve the 
following problem: Suppose that we are given an oper- 
ator O acting on an N dimensional Hilbert space and a 
quantum state p. By this we mean that someone supplies 
us with many copies of a quantum system prepared in the 
same state p and defines for us the operator O by spec- 
ifying its expansion in a basis of the space of operators. 
Our task is to compute the expectation value of O in 
the state p. We will show that it is possible to construct 
a programmable circuit that evaluates such expectation 
value by measuring the polarization of a single qubit. 
The inputs of such circuit are a data register, a program 
register and an auxiliary qubit. The circuit evaluates the 



expectation value of an operator O (specified by the pro- 
gram register) in the quantum state p of the data register. 
The expectation value Tr(pO) is obtained by performing 
a measurement of the polarization of the auxiliary qubit. 
We will describe how to construct these circuits and ex- 
hibit an interesting example: a programmable array to 
efficiently solve a class of quantum decision problems con- 
cerning properties of quantum phase space distributions. 

The quantum gate arrays discussed in this letter are 
designed using the "scattering circuit" shown in Figure 1 
as a simple primitive. In such circuit a system, initially in 
the state p, is brought in contact with an ancillary qubit 
prepared in the state |0) . This ancilla acts as a probe par- 
ticle in a scattering experiment. The algorithm consist 
of the following steps: i) Apply a Hadamard transform 
H to the ancillary qubit. Since H\0) = (|0) + |1))/V2, 
H\l) = (|0) - |1))/V2, the new state of the qubit is 
(1°) + |l))/\/2. ii) Apply a "controlled-^" operator, 
which does nothing if the state of the ancilla is |0) and 
applies the unitary operator A to the system if the an- 
cilla is in state 1 1) , iii) Apply another Hadamard gate to 
the ancilla and perform measurements of its spin polar- 
izations along the z- and y-axes. Given sufficiently many 
independent instances of the experiment, the measure- 
ments yield the expectation values (<j z ) and (a y ) of the 
Pauli spin operators a z and a y . This algorithm has the 
following remarkable property: 



(a z ) = Re[Tr(Ap) ], (a y ) = lm[Tr{Ap)[ 



(1) 



Different versions of this circuit play an important role 
in many quantum algorithms [a, 13, lj- 111 ■ In particular, 
the scattering circuit was recently used as a basic tool 
to interpret tomography and spectroscopy as two dual 
forms of the same quantum computation |ll| . 

For our purpose it is useful to review how to use this 
scattering circuit as a primitive to design a tomographer 
(i.e., a device that after a number of experiments deter- 
mines the quantum state of the system). As a conse- 
quence of equations JJ| we see that every time we run 
the algorithm for a known operator A, we extract in- 
formation about the state p. Doing so for a complete 
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FIG. 1: The scattering circuit that can be used to evaluate 
real and imaginary parts of the expectation value Tr(pA) for 
a unitary operator A. H denotes a Hadamard transform. 
The "controlled- A" operation is such that (ctrl— A) \q) = 
\q)A«\9). 



basis of operators {A(a)} one gets complete information 
and determines the full density matrix. Different tomo- 
graphic schemes are characterized by the basis of oper- 
ators A(a) they use. Of course, completely determining 
the quantum state requires an exponential amount of re- 
sources. In fact, if the dimensionality of the Hilbert space 
of the system is N then the complete determination of the 
quantum state involves running the scattering circuit for 
a complete basis of iV 2 operators A(a). However, evalu- 
ating any coefficient of the decomposition of p in a given 
basis can be done efficiently provided that the operators 
A(a) can be implemented by efficient networks. A con- 
venient basis set is defined as (see, for example, [TlTll^l: 



A(a) = A(q,p) = U q RV~ p exp(iirpq/N). 



(2) 



Here, both q and p are integers between and N — 1, 
U is a cyclic shift operator in the computational ba- 
sis (U\n) = \n + 1)), V is the cyclic shift operator in 
the basis related to the computational one via the dis- 
crete Fourier transform, and R is the reflection operator 
(R\n) = \N — n)). It is straightforward to show that 
the operators A(q,p) are hermitian, unitary and form a 
complete orthonormal basis of the space of operators sat- 
isfying 

Tr[A(a)A(a')} = N 6 N (q' - q)5 N {p' - p), (3) 

where <Jjv(x) is the periodic Kronccker delta function that 
is equal to one if x = (modulo N) and vanishes oth- 
erwise (the above operators form a "quorum", as de- 
fined in [Hj). With this choice for A(q,p) the scatter- 
ing circuit directly evaluates the discrete Wigner function 

We will now show how to design a programmable gate 
array to evaluate the expectation value of any operator O. 
We will assume that we know how to expand O in a ba- 
sis such as the one used above: O = ^2 q p o(q,p)A(q,p). 
As the operators A{q 1 p) are not only unitary but also 
hermitian, the real and imaginary parts of the complex 



coefficients o(q,p) define the expansion of the hermitian 
and anti- hermitian pieces of O in the basis A(q,p). The 
expectation value of these two pieces can be evaluated 
separately using the procedure described below (the re- 
sults can then be combined to get the expectation of O). 
So, in what follows we will assume that the operator at 
hand is hermitian and that the coefficients of its expan- 
sion in the basis A(q,p) are real numbers. To introduce 
our method, it is convenient to notice first that the eval- 
uation of the expectation value of the operators A(q,p) 
can be done using a programmable circuit that is inde- 
pendent of q and p. Such circuit is illustrated in Figure 
2. This is an application of the scattering circuit shown 
in Figure 1 with two program registers used to encode 
the value of q and p. When the quantum state of the 
program is \^f)p = \q)\p) the circuit evaluates the expec- 
tation value of A(q,p). This is accomplished by letting 
the program registers to act as controls of the operators 
U and V (which, as mentioned, generate cyclic shifts ei- 
ther in the computational or in the conjugated basis). 
Thus, the action of the circuit is such that when the aux- 
iliary qubit is in state \a) {a = 0, 1) and the state of the 
program is \q}\p), the operator A a (q,p) is applied to the 
system register. The network in Figure 2 is efficient since 
it can be built using a number of elementary gates which 
scales polynomially with log(iV) [l4j |. 

The circuit has an obvious property: Different states 
\q)\p) are used to program the evaluation of the expec- 
tation value of orthogonal operators A(q,p). It is clear 
that by restricting to such program states one has no 
real advantage with respect to the case in which q and p 
are stored as classical information. However, we can use 
more general program states: If the program register is in 
the state \^)p = J2 q p c (.1^P)\l)\p) then the same circuit 
evaluates the expectation value of a linear combination 
of the operators A{q 1 p) since the final polarization is: 



(a z )=Tr[p^2\c(q,p)\ 2 A(q,p) 



(4) 



Equation J3J shows that this algorithm can be used to 
evaluate the expectation value of any operator that can 
be written as a convex sum of the basis set A(q,p). This 
is not general enough since the expansion of a hermi- 
tian operator can include negative coefficients. For this 
purpose the method can be extended as follows: The 
most general hermitian operator can be expressed in the 
basis A(q,p) as O = J2 q , p c 2 (q,P) Gxp(i7T(j>(q,p))A(q,p), 
where c(q,p) is a real number and (j>(q,p) is either or 1 
(4>(q,p) simply stores the information about the sign of 
each coefficient). We will assume that J2 q p c 2 {q,p) = 1 
(if this is not the case we can always renormalize the co- 
efficients). For operators with some negative coefficients 
we can use a third register consisting of a single qubit 
to store 4>(q,p). The circuit evaluating the expectation 
value of O is shown in Figure 3. The first two program 
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FIG. 2: Programmable gate array evaluating Tv(pA(q,p)) 
from the polarization of the first qubit. The program state 
is \^)p = \q)\p}- All "controlled-O" operators act as: 
(ctrl-0)|n)|*) = \n)O n \^). U K (V K ) are cyclic shift op- 
erators in the computational (conjugated) basis of a K di- 
mensional space. A subscript in an operator denotes the di- 
mensionality of the space in which it acts 
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FIG. 3: Programmable gate array evaluating Tr(pO) from 
the polarization of the first qubit. The program state is 

= I2 q ,p c (<l>P)\<l}\P}\®(<liP)) where c(q,p) and 
define the polar decomposition of the coefficients o(q,p) = 
Tr(OA(q,p))/N. 



registers store q and p and are used exactly in the same 
way as above. The third one, storing cj>(q,p), is acted 
upon with a a z operator, introducing the required phase 
exp(iir(j)(q,p)). Then, if the state of the program register 

is I*)p = Eg,p c (9 J -P)l'?)b)l ( / , (9 1 -P)>- the final polariza- 
tion measurement turns out to be 



Tr 



V Q,P 



{q,p)e l7T ^ q '^A(q,p) = Tr(pO) . 



(5) 

Summarizing, we showed that the measurement of the 
expectation value of any operator O can be done using 
a programmable gate array. The hardware architecture 
is associated with the particular choice of basis A(q,p), 
which is just a matter of convenience, and is indepen- 
dent of O. The software used to program the array is 
obviously determined by the choice of hardware. The 



expectation values of the hermitian and anti-hermitian 
parts of the operator O are computed separately using 
a method that requires knowledge of the expansion of 
these operators in the basis A(q,p). If the coefficients in 
the expansion are written as o(q,p) = Tr (OA(q,p))/N — 
c 2 (q,p) exp(i7T0(g,p)), the program state that needs to 
be prepared is \9) P = E g , P c fop) k)\p)\Hq,P)) (where 
4>(q,p) — or 1). Thus, the coefficients of the expansion 
of O in the basis A(q,p) define the program state \^f)p 
required to measure its expectation value. It is clear 
that in most cases the method will not be efficient. For 
example, both the task of defining the operator by speci- 
fying the coefficients o{q 1 p) as well as the preparation of 
the program state \^)p are likely to be inefficient. The 
existence of efficient networks to implement "controlled- 
A(q,p)" operations is a less stringent conditions that is 
fulfilled by the basis defined in @ 14] . Having said this, 
it is worth noticing that there are sets of problems that 
can be efficiently solved using this method. We will now 
describe one such example. 

We will show that the circuit of Figure 2 can be easily 
adapted to evaluate the sum of values of the Wigner func- 
tion over various phase space domains. The program reg- 
ister is used to define the domain over which the Wigner 
function is averaged. If the domain is a line, the algo- 
rithm just evaluates the probability for the occurrences of 
the results of the measurement of a family of observables 
(see below). However, for more general domains (such as 
line segments, parallelograms, etc) the circuit evaluates 
properties that characterize a quantum state that cannot 
be simply casted in terms of probabilities. In this sense 
the circuit is as a programmable tomographer measur- 
ing various features of phase space distributions. Before 
going into more details let us briefly review some proper- 
ties of Wigner functions. Discrete Wigner functions can 
be used to represent the quantum state of a system in 
phase space 0, 0, 0, Q|. For a system with an TV 
dimensional space of states such function is defined on a 
lattice of 2N x 2N points (q,p) where both q and p take 
values between and 2N — 1 (only N x N oi these values 
are independent) . At each phase space point the Wigner 
function is defined in terms of the operators A(q,p) given 
in © as 



W(q,p) = ^Tr(A(q,p)p). 



(6) 



As mentioned above, the measurement of this function 
can be done by using the scattering circuit 01 or the pro- 
grammable gate array of Figure 2. The program register 
encodes the value of q and p. In general, if the program 
state is \^)p = p c{q,p)\q)\p) , the final polarization 
measurement is (a z ) — 2N^ qp \c(q,p)\ 2 W{q,p). Thus, 
the program defines the region over which we average 
the value of the Wigner function. In general, prepar- 
ing the program state associated with a general phase 
space region can be complicated. However, there are 
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simple procedures to prepare the program states corre- 
sponding to general lines, segments and parallelograms. 
Let us begin with the simplest case: For the program 
state \^}p = J2 q \q) \po)/y2N, it is clear that the final 
polarization measurement reveals the sum of values of 
the Wigner function along the horizontal line defined as 
p = pq: i.e., (a z ) — J2 q W(q,po). To consider more gen- 
eral lines one can notice that the horizontal line p = c 
can be mapped into one satisfying the equation p — bq = c 
(mod 27V, with b and c integers between and 27V) by 
applying a linear area preserving map. Unitary opera- 
tors corresponding to quantizations of such maps (known 
as "quantum cat maps") act classically in phase space. 
Thus, Wigner function flows from one point to another 
according to the classical transformation. Using this fea- 
ture, one can compute the sum of the Wigner function 
along any tilted line as follows: One first transforms the 
state with the appropriate cat map and later evaluates 
the sum of the Wigner function along a vertical or hor- 
izontal line. Moreover, the method can be made fully 
programmable by adding extra registers to store the in- 
tegers b and c parametrizing any cat map. Evaluating 
sums of Wigner functions over phase space lines is par- 
ticularly interesting because of a crucial property of such 
functions: Thus, adding W(q,p) along the line ap—bq = c 
one obtains the probability to detect the eigenstate of the 
translation operator T(b,a) — U a V b exp(inab/N) with 
eigenvalue exp(z7rc/7V). As a consequence, in this case 
the programmable gate array evaluates probabilities for 
the possible results of a set of measurements. Lines are 
a special case since more general phase space domains 
cannot be associated with projection operators. How- 
ever, it is clear that the method described above can be 
applied to efficiently prepare program states for other 
phase space domains such as general (tilted) parallelo- 
grams: One first trivially prepare the program state for 
a parallelogram limited by vertical and horizontal seg- 
ments and later tilt it applying the strategy based on the 
use of cat maps. Other simple phase space regions can 
also be programmed using variations of this method. It 
is also interesting to notice that using variations of the 
circuit shown in Figure 3 we can also subtract values 
of the Wigner function in different phase space regions 
(which could be useful if one is interested in comparing 
their values). 

In this letter we established the existence of a gate ar- 
ray that can be programmed to evaluate the expectation 
value of any operator acting on an TV dimensional Hilbcrt 
space. The expectation value is obtained by measuring 
the polarization of a single auxiliary qubit. As an exam- 
ple, we showed how to program the evaluation of sums of 
values of the discrete Wigner function over various sim- 
ple phase space domains. It is important to mention that 
our method is only efficient to determine if the sum of 
the Wigner function in a phase space domain (with up 
to o(7V) points) is greater than a fixed, TV-independent, 



threshold (since this does not require exponential preci- 
sion) . This is a "quantum decision problem" whose input 
data (encoded in the system's state p) is inherently quan- 
tum. Due to the nature of the input data, this problem 
cannot even be formulated on a classical computer. Inter- 
est on problems with quantum input data have recently 
increased, partly due to their significance in connection 
with the potential detection of entanglement 0, 0, 
as well as their relation with tomographic problems like 
the one described here. The extension of some of the 
above results to continuous variables is still under investi- 
gation [J3 ■ After completing this work we became aware 
of the related approach to the construction of quantum 
universal detectors presented in (lif. 
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